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by a structure subjected to echographic examination, 
comprising the steps of: a) transmitting an ultrasonic ex- 
citation signal into a portion of said structure under ex- 
amination; b) receiving a radio frequency response sig- 



nal from said structure; c) applying a time-frequency 
transform to said radio frequency response signal, di- 
viding the radio frequency response signal into a plural- 
ity of frequency bands; d) calculating a local spectral pa- 
rameter from the values of the time-frequency trans- 
form. 
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Description 
Technical field 

5 [0001] The present invention relates to a method for processing echographic signals, particularly in the field of non- 
invasive clinical examination based on the ultrasonic method. The invention also relates to an ultrasonic device for 
applying a method for processing echographic signals. 

Prior art 

10 

[0002] Echography is a non-invasive method of investigation for examining internal organs of the human and animal 
body, and also for the non-invasive examination of other types of structures, including non-living and inorganic struc- 
tures. As is known, it is based on the generation of a set of ultrasonic signals (forming an acoustic beam) which are 
directed into the patient's body and which excite an ultrasonic response signal, caused by back reflection and/or dif- 
15 fusion of sound waves by the diffusers and/or reflectors present within the region under observation. The probe for 
transmitting and receiving the echographic signal sends the response signal to a receiving and processing system, 
after which the signal is converted into a two-dimensional ultrasonic image. 

[0003] Commercial ultrasonic scanners supply images relating to the amplitude of the echographic signal back- 
propagated by the various tissues or interfaces struck by the acoustic beam. 

20 [0004] The various tissues or the various interfaces encountered are detected on the basis of their echogenicity 
determined by their acoustic impedance which in turn is the product of the density and the acoustic velocity. 
[0005] The purpose of conventional ultrasonic scans is to reconstruct and maintain the spatial ordering of the tissue 
under examination, by making use of amplitudes determined in the time domain, which in turn yields the spatial location 
in depth. The signal treatment which is carried out has the aim of maximizing the signal to noise ratio and holding the 

25 sensitivity of the examination constant while the depth varies. 

[0006] When the acoustic wave is propagated in a medium and encounters a non-homogeneous region, some of 
the energy is back-propagated towards the transducer and some is propagated into the tissue. 
[0007] The back-propagated acoustic signal consists of the temporal superimposition of the contributions from each 
reflector and/or diffuser; the result is a signal having an extremely complex shape, due to the random nature of the 

30 distribution and the shapes and mechanical properties of the diffusers and the reflectors. The acoustic signal back- 
propagated by the medium under examination is converted by the transducer into an electrical signal, called the radio 
frequency signal. 

[0008] Fig. 1 shows a portion of a radio frequency signal forming part of an ultrasonic trace relating to biological 
tissue. More particularly, Fig. 1 shows a portion of a track of the return echographic signal from a portion of an organ 
35 subjected to echographic examination. The horizontal axis shows the response time, in other words the depth of the 
location of the diffuser or reflector which generates the return signal captured by the transducer. The vertical axis shows 
the amplitude of the signal. The ultrasonic image of the organ under examination is reconstructed from a set of tracks 
of this type, which form a frame of signals. 

[0009] A first very simple processing operation carried out on the radio frequency signal is the extraction of the signal 
40 amplitude by means of an envelope extraction operation. In Fig. 2, the letter I indicates the result of the envelope 

operation carried out on the radio frequency signal shown in Fig. 1 . It will be noted that the envelope signal follows the 

positive crests of the radio frequency signal and therefore provides an estimate of the signal amplitude. 

[0010] The envelope operation is currently used in all ultrasonic imaging techniques and forms the basis of the A- 

mode presentation of the individual trace (where A stands for amplitude) and B-mode presentation (where B stands 
45 for brightness) for the resulting images, each individual A-mode trace being encoded on a gray or color scale. 

[0011] In ultrasonic scans of the present generation, the envelope of the echographic signal is used, in other words 

the detected signal is sent to the A-mode, B-mode and M-mode traces. 

[0012] With this solution, there is a loss of the information contained in the radio frequency signal, in other words the 
information in the form generated by the receiving transducer, consisting of a signal having a carrier frequency approx- 
50 imately equal to that of the center frequency of the band of said receiving transducer, with a frequency modulation 
(asymmetric spectrum) with phase steps (frequency peaks) and possibly with a morphology due to non-linear distortion 
when the medium shows non-linear effects (generation of harmonic spectra). 

[0013] Fig. 3 shows the envelope signal of the radio frequency trace of Fig. 1 , in a conventional A-mode represen- 
tation. It will be noted that information on the frequency and phase of the original signal has been lost irretrievably. 
55 [0014] In particular, if the radio frequency is eliminated by an envelope extraction operation, there is a loss of: 

• information obtained from non-linear interactions between ultrasound and biological tissues (local acoustic imped- 
ance which depends on the density of the medium, which in turn depends on the intensity of the incident ultrasonic 
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wave, re-emission on the harmonics and subharmonics by bubbles present in the contrast medium, etc.) 
• information obtained from linear interactions, namely the modifications of the signal spectrum due to selective 
absorption, especially where this is depth-dependent, and the scalloping effect (indentation of the shape of the 
amplitude spectrum present in signals which have been backscattered by discontinuities, including those of the 
5 point type, distributed periodically in space) 

[001 5] The conventional echography systems used for medical diagnosis are therefore based on the observation of 
an A-mode trace or a B-mode image, which show only the amplitudes of the echo signal relating to the various tissue 
structures encountered. 

10 [0016] Consequently there are situations in which the conventional methods are found to be limited since, owing to 
low sensitivity or low specificity, it is not possible to obtain a reliable diagnosis of pathology. However a radio frequency 
signal contains not only the information relating to amplitude, but also the frequency and phase information (related to 
the time delays with which the echoes due to individual diffusers reach the transducer). 

[0017] The latest signal processing methods tend to extract all the parameters characterizing the signal, and in 
15 particular what is known as the spectral content. Spectral analysis applied directly to the radio frequency signal can 
extract information from all the characteristics of the echo signal as it is returned from the transducer. 
[0018] A spectral analysis algorithm for "reading" the signal in the frequency domain makes it possible to provide 
characteristic values for the differentiation of the various types of signal and consequently of the biological structures 
with which they interact. 

20 [0019] Diagram (A) in Fig. 4 again shows the radio frequency signal of Fig. 1 , while diagrams (B) and (C) show the 
spectral content of the signal for two portions of the trace of Fig. 4A, to show how the structural differences in the 
tissues result in different types of ultrasonic impulses and consequently distinctive frequency variations in the corre- 
sponding characteristics. The spectral contents of the two portions of signal shown in Figs. 4B and 4C extend from 0 
to 12 MHz. It should be noted that both spectra have approximately the same central frequency (the frequency at which 

25 the peak is found) in the vicinity of 4 MHz (the signal relates to scanning with an ultrasonic probe having a nominal 
frequency of 5 MHz), but have different frequency variations, since the first graph shows contributions (above 5 MHz) 
which are present to a lesser extent in the other portion of the signal. 

[0020] The above facts demonstrate that spectral analysis of the radio frequency signal is a powerful tool for char- 
acterizing tissues by means of the corresponding ultrasonic signal. 
30 [0021] Recently, the development of new signal processing algorithms has made it possible to obtain further infor- 
mation on tissues, the trend being to investigate the interaction with the ultrasonic beam in terms of spectral variations 
induced in the reflected signal. In fact, it has been found that pathological structures can be identified according to the 
differences in their selective frequency filtering action. 

[0022] Examples of attempts to type tissues by spectral analysis of ultrasonic signals are described in E.J. Feleppa 
35 et al, Typing of Prostate Tissue by Utrasonic Spectrum Analysis, in IEEE Transactions of Ultrasonics, Ferroelectrics 
and Frequency Control, Vol. 43, no. 4, July 1996, p. 609 ff.; E.J. Feleppa et al., Diagnostic Spectrum Analysis in 
Ophthalmology: A Physical Perspective, in Ultrasound in Med. & Biol. Vol. 12, no. 8, 1986, p. 623 ff.; G. Schmitz, et 
al., Tissue-Characterization of the Prostate Using Radio Frequency Ultrasonic Signals, in I EE Transactions on Ultra- 
sonics, Ferroelectrics and Frequency Control, Vol. 46, no. 1, January 1999, p. 126 ff.; F.L. Lizzi, Utrasonic Scatterer- 
40 Property Images of the Eye and Prostate, in 1997 IEEE Ultrasonics Symposium, p. 1109 ff.; E.J. Feleppa et al, 3-D 
Tissue Typing of Prostate Tissue Based on Spectral Parameters, in 1 995 Ultrasonics Symposium, p. 1 1 71 ff . ; G. Schmitz 
et al, Tissue Characterization of the Prostate Using Kohonen-Maps, in 1 994 Ultrasonics Symposium, p. 1487 ff.; A.L. 
Huynen et al., Analysis of Ultrasonographic Prostate Images for the Detection of Prostatic Carcinoma: The automated 
Urologic Diagnostic Expert System, in Ultrasound in Med. & Biol., Vol. 20, no.1 ., 1994, p. 1 ff. 
45 [0023] However, the methods developed up to the present time have not proved to be effective, either because they 
are unsuitable for typing tissue in a reliable way, or because of the excessive computing capacity required for the 
algorithms used, which cannot be implemented in a way that permits the real time display of the results, which would 
make it possible to envisage clinical applications where the speed of data output is a basic requirement. 

50 Objects and brief description of the invention 

[0024] The object of the present invention is to provide an effective and reliable method for typing a tissue by the 
spectral analysis of an echographic signal. More generally, the object of the present invention is to provide a method 
for characterizing a generic structure which may be an organ of a human or animal body, but which may also be an 
55 inanimate composite structure. This is because the method according to the present invention can be applied both to 
ultrasonic analysis in medicine and diagnosis and to more purely technological fields, for carrying out non-invasive 
examinations to determine the characteristics of a structure, for example in order to check its internal integrity. In the 
present description and in the attached claims, therefore, the term "structure under examination" denotes, in a general 
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way, either a tissue or a living organ, or an inanimate structure. 

[0025] This and further objects and advantages, which the following text will make clear to those skilled in the art, 
are achieved essentially by means of a method of spectral analysis of an ultrasonic signal returned by a structure 
subjected to echographic examination , comprising the steps of: 

5 

a) transmitting an ultrasonic excitation signal into a portion of said structure subjected to examination; 

b) receiving a response signal from said structure; 

c) applying a time-frequency transform to said response signal, thus dividing the response signal into a plurality 
of frequency bands; 

10 d) calculating a local spectral parameter from the values of the time-frequency transform, by statistical investigation 

of the values of the time-frequency transform. 

[0026] The time-frequency transform is an essential requirement for the local extraction of the information contained 
in the frame of the back-propagated echographic signal. 

15 [0027] The method of the present invention is essentially based on the determination, for a number of samples of 
the input frame, of a parameter which takes into account the variations in the characteristics of the signal spectrum. 
The values of this parameter are arranged in a frame which, in a subsequent display stage, generates a color scale, 
which can be superimposed on the B-mode image of the ultrasonic scan, and which provides an indication of the 
variation of the spectral characteristics of the signal. Since it has been found that pathological tissue structures, for 

20 example tumor cells, cause an alteration of the spectral characteristics of the ultrasonic return signal, this is reflected 
in an alteration of the value of the local spectral parameter. An indication of the location of any neoplasms of the 
examined organ can be obtained and superimposed on the conventional B-mode image, by using a color scale to 
display the areas of the image in which the spectral parameter statistically takes values lying within a reference class 
which, experimentally, can be associated with neoplastic phenomena. 

25 [0028] The invention, which is based on the considerations stated above regarding the loss of information due to the 
operation of extracting the envelope of the radio frequency signal, essentially makes it possible to detect the local 
information for the structure under examination by the extraction of the local spectrum either within a small time slot, 
or within a small range of distances around each point of the echogram. 

[0029] In practice, the method is based on the principle of characterizing the echographic signal by decomposing it 
30 into frequency bands by means of a time-frequency transform, using a suitable decomposition algorithm, and proceed- 
ing to analyze the variation of the coefficients of the time-frequency transform at each instant of time in the various 
frequency bands. 

[0030] Theoretically, it is possible to use any time-frequency transform which can encode the suitably sampled and 
digitized echographic signal into bands or sub-bands. The frequency band encoding makes it possible to identify the 
35 variation of the coefficients in the various bands, and thus essentially to identify the frequency variations of the signal 
for individual instants of time. The local spectral parameter is extracted from these variations. This operation can be 
carried out by using the convolution of the backscattered signal with a special waveform (wavelet) which has a mor- 
phology similar to that of the transmitted signal, but has a continuously variable carrier frequency, on the signal whose 
spectrum is to be determined. 

40 [0031] The point spectrum which is obtained has a morphology which depends on the linear and non-linear interaction 
effects mentioned above; these effects are due to the mechanical characteristics and microstructural organization (in 
the dimensional scale around the values of the wavelength of the interacting ultrasounds) of the small elementary 
volume (point) in question. For faster calculation, the discrete form of wavelet transform is preferably used. It has also 
been found that appreciable results are obtained by using the discrete wavelet packet transform. 

45 [0032] The variation of the local spectrum is approximated by a polynomial function whose coefficients are organized 
in classes of values. 

[0033] The classes are made significant by examining the values of the coefficients derived from spatial means in 
intervals along the line of sight around the point in question and on corresponding portions of a certain number of lines 
of sight adjacent to both sides of the line of sight on which the point in question is located. 

50 [0034] This is because the decomposition algorithm based on the discrete wavelet packet transform (DWPT) makes 
it possible to have high frequency resolutions in each region of the signal band, and also to analyze the time-frequency 
plane with a different resolution. For example, after fixing the frequency region in which higher resolution is desired, 
the corresponding sequence of coefficients is decomposed at a higher level, while the others can be decomposed at 
a lower level, with advantages in terms of temporal resolution. 

55 [0035] In a practical embodiment of the invention, the local spectral parameter is calculated on the basis of a frame 
of the energy content of the signal, found for each band into which the digitized signal has been sampled. 
[0036] In practice, the method according to the present invention can comprise the steps of: 
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acquiring an input frame, consisting of a plurality of tracks, each track consisting of the radio frequency signal 
received from the transducer after the generation of an excitation pulse; 

sampling and digitizing each track of said frame, thus generating a digitized frame consisting of a plurality of tracks, 
each containing a set of digitized samples; 
5 - applying said time-frequency transform to the digitized frame, thus encoding said digitized frame into a plurality of 
frequency bands, generating one frame of the transform in each band, 

calculating the energy content of the signal for each track of the frame of the transform in each band, thus defining 
an intermediate energy content frame, said energy content frame having a set of samples representing the energy 
content of the input signal for each track and in each band; 
10 - interpolating for each track of the energy content frame, with an interpolating polynomial above the second degree, 
the samples of all the bands for each instant of time, one of the coefficients of said interpolating polynomial rep- 
resenting said local spectral parameter, and generating a frame of the spectral parameters containing the values 
of said coefficient for each coefficient of the energy content frame. 

15 [0037] More generally, the method can include encoding the digitized signal into frequency bands, and consequently: 

selecting at least two frequency bands into which said response signal has been divided; 
extracting a local spectral parameter from the coefficients of said at least two frequency bands. 

20 [0038] Essentially, the method according to the present invention is based on the concept of decomposing the signal 
into frequency bands and then extracting, from the variation of the coefficients in the individual bands, local statistical 
information related to the characteristics of the tissues or more generally of the structure which has generated the 
echography response signal. In practice, the following steps may be provided: 

25 - acquiring an input frame, consisting of a plurality of tracks, each track consisting of the response signal corre- 
sponding to an excitation signal; 

sampling and digitizing each track of said frame, thus generating a digitized frame consisting of a plurality of tracks, 
each containing a set of digitized samples; 

applying said time-frequency transform to the digitized frame, thus encoding said digitized frame into a plurality of 
30 frequency bands, generating one frame of the transform in each band; 

generating from the coefficients of at least two of said frames of the transform a frame of local spectral parameters, 
in which each coefficient is obtained by processing the coefficients of said at least two frames of the transform. 

[0039] The local spectral parameters can be found by various forms of processing, according to the type of local 
35 information which is to be extracted from the signal encoded into frequency bands. For example, it is possible to 
determine a spectral parameter related to the variation of the phase of the back-propagated signal. 
[0040] Alternatively, it is possible to move a window over each of the two frames of the transform, and to determine, 
for each step, the maximum value taken by the coefficients in said window. The maximum values extracted for each 
step in the two frames of the transform can then be combined with each other to obtain the frame of the local spectral 
40 parameters. Each frame can then be subjected to statistical investigation, for example in order to determine the pop- 
ulation, in other words the statistical distribution of the parameters in one or more reference classes. 
[0041] By determining the maximum values taken by the coefficients for each position of the window as it moves 
over the frames of the transform, an intermediate frame is obtained for each frame of the transform. 
[0042] When the processing which is carried out is based on the energy content of the signal, the intermediate frame 
45 is represented, for each frame of the transform, by the corresponding energy content frame. 

[0043] Further advantageous characteristics and embodiments of the method according to the invention are indicated 
in the attached claims and are described in greater detail below with reference to a practical embodiment of the method. 
[0044] The invention also relates to an echography device comprising an echography probe and means of acquiring 
and processing a return signal from an organ subjected to echographic examination, characterized in that said acqui- 
re sition and processing means are programmed to implement a method as specified above. 

Brief description of the drawings 

[0045] The invention will be more clearly understood from the description and the attached drawing, which shows a 
55 practical and non-restrictive embodiment of the invention. More particularly, 

Figs. 1 to 4, described above, show an echography trace and a representation of the spectral content in two portions 
of said trace; 
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Fig. 5 shows a simplified block diagram of an echographic apparatus; 

Fig. 6A shows an echographic excitation signal sent by the transducer into the organ to be examined; 
Figs. 6B to 6D show, respectively, the RF response signal , the Daubechies 1 6 wavelet used forf iltering the digitized 
RF signal, and the coefficients of the DWPT obtained by decomposition into eight bands; 
5 Figs. 7A-7G show schematically the sequence of processes carried out on the echographic signal by the method 

according to the present invention; 

Fig. 8 is a block diagram summarizing the sequence of processes; 

Fig. 9 shows a filter circuit for the decomposition of the signal by the discrete wavelet packet transform (DWPT); 

Fig. 10 shows an implementation of the decomposition algorithm with a single pair of filters; 
10 Fig. 11 shows an implementation of the same decomposition algorithm, using parallelized filters; 

Fig. 12 shows the band coefficients obtained by application of the DWP transform for a track of the input RF signal; 

Fig. 13 shows schematically the process for calculating the signal energy; 

Fig. 14 is a block diagram of the module for calculating the signal energy; 

Fig. 1 5 shows the flow chart for the calculation of the signal energy; 
15 Fig. 16 shows schematically the interpolation operation for determining the local spectral parameter for a given 

track of the signal; 

Fig. 1 7 shows the block diagram for calculating the local spectral parameter; 
Fig. 1 8 shows the flow chart for the calculation of the local spectral parameter; 

Fig. 1 9 shows the block diagram for classifying the coefficients relating to the spectral parameter and for generating 
20 a frame containing the frequency of the distribution of the spectral parameter in the classes examined; 

Fig. 20 shows the flow chart corresponding to the block diagram of Fig. 1 9; 

Fig. 21 shows the block diagram for calculating the mean of the coefficients of the frame containing the frequency 
of the distribution of the spectral parameter; 

Fig. 22 shows the flow chart corresponding to the block diagram of Fig. 21 ; 
25 Fig. 23 shows schematically a different implementation of the method according to the invention. 

Detailed description of the preferred embodiment of the invention 

General description of the apparatus 

30 

[0046] Fig. 5 shows a highly schematic block diagram of an echographic apparatus. The number 1 indicates a trans- 
ducer of the scanning head, indicated in a general way by 2, which emits the echographic excitation signal and receives 
the echographic response signal. The analog output 3 of the scanning head is connected to an acquisition card 5. This 
generates a sampled and digitized radio frequency signal, which forms the input for a bank of integrated FIR filters, 

35 indicated in a general way by 7. The bank of filters 7 can consist of any commercially available integrated device, for 
example the 431 68 circuit produced by Harris, the GF 1 91 circuit produced by Gennum, the GEC 1 6256 circuit produced 
by Plessey, or the GC2011 circuit produced by Graychip. The FIR filter 7 is associated with a digital signal processor 
(DSP) indicated by 9, which acts as a controller for the whole processing unit. The microcontroller can be, for example, 
a TMS 32031 DSP made by Texas (USA) or other equivalent. 

40 [0047] The number 11 indicates a store associated with the filter 7 and with the controller 9. The store 1 1 is connected 
to an analog/digital converter 13 and to the digital stages of the scanning head 2. 

[0048] The circuit described here is known and does not require more detailed description. Essentially, the apparatus 
operates as follows. The transducer 1 emits a series of echographic excitation signals having, for example, the form 
shown in Fig. 6A, into the body containing the organ to be examined, for example the prostate, the liver or other part. 
45 The tissue struck by the ultrasonic waves of the signal emitted by the transducer 1 returns a scattered or reflected 
signal which has the form of the signal in Fig. 1 (radio frequency signal). 

[0049] The radio frequency signal is taken from the receiving circuit of the echograph immediately after the time gain 
compensation (TGC) unit and is then digitized at a suitable sampling frequency, for example 40 MHz, and with a 
resolution of 12 bits (or other suitable resolution), by the acquisition card 5. In practice, acquisition is carried out for 

50 individual tracks until a frame or a portion of a frame is complete. The signal samples acquired for each individual track 
are stored in a memory bank of the acquisition card until a complete frame or a portion of a frame has been acquired; 
it is only at this point that the stored data are sent to the processor card comprising the DSP 9 and the FIR filter 7 which 
execute the operations which are described in detail below, according to the method of the present invention. While 
the data is being transferred to the processor card, the acquisition card starts to acquire the signal samples of the 

55 tracks for the next frame. The processing time of the processor card 12 is short enough to allow the processing to be 
completed within the acquisition time of the sampled signal for the next frame. Thus real time operation is achieved. 
[0050] The processorcard 12carries out the processing of thesignal, particularly the time-frequency transform using 
the integrated filters 7 and the other operations which are described in greater detail below, to obtain the characterization 
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of the signal spectrum. 

Brief description of the signal processing procedure 

5 [0051] Figs. 7A-7E and 8 summarize, in the form of graphs and a summary flow chart, the operations carried out on 
the input signal to obtain the representation of the frame containing the information on the local spectral parameter 
determined by the method according to the present invention. 

[0052] The processing of the sampled and digitized signal acquired by the card 5 commences with the application 
of a time-frequency transform to encode the signal into sub-bands. Essentially, for each track of the acquired, sampled 

10 and digitized frame, the signal is decomposed into a series of frequency bands or sub-bands by the application of a 
wavelet transform, and in particular by means of the discrete wavelet packet decomposition. The filtering is performed 
by the bank of integrated filters 7. The number of bands or sub-bands into which the digitized signal is decomposed 
depends on the number of consecutive filtering operations used. In practice, a decoding into eight bands or sub-bands, 
indicated below as band 0, band 1 , ... band 7, is sufficient for the purposes of the method according to the present 

15 invention. 

[0053] Fig. 7 A shows schematically a frame acquired by the scanning head of the echograph. The frame, referred 
to below as the "input frame", consists of a plurality of n tracks. Each track consists of an analog signal represented 
by the wavefront reflected or back-propagated by the tissues in which the ultrasonic excitation signal is propagated. 
After the sampling and digitization, we have a frame (Fig. 7B) which also consists of n tracks, each of which contains 

20 a series of digitized signal samples. This frame will be referred to as the "digitized frame". 

[0054] By applying the time-frequency transform to the digitized frame, in other words by filtering the sampled and 
digitized signal consisting of the set of coefficients forming the digitized frame (shown schematically in Fig. 7B), we 
obtain a distribution in bands or sub-bands (eight in the present example): in practice, for each track from 1 to n there 
is, for each instant of time, in other words for each sample in the time direction (in other words in the direction of 

25 propagation of the signal), a coefficient in each of the eight bands shown schematically in Fig. 7C and indicated as 
band 0, band 1 , ... band 7. The coefficients of the time-frequency transform applied, in this case the discrete wavelet 
packet transform (DWPT), are shown in the various bands. The frames containing the coefficients of the time-frequency 
transform are referred to below as the "transform frames". 

[0055] In a subsequent step, the coefficients of the transform frame, in the various frequency bands, are processed 
30 to extract a distribution of the energy content of the signal. This operation is carried out on each track of the frame and 
in each band. It essentially consists in the construction of a frame, referred to below as the "energy content frame", 
whose coefficients represent the energy content of the signal along each track, in other words in the direction of the 
depth of propagation of the signal. Fig. 7D shows schematically the coefficients representing the energy content of the 
signal, which are divided into the eight frequency bands and which form the energy content frame for each band. 
35 [0056] The next operation consists in the extraction of a local spectral parameter of the signal for each point of the 
energy content frame (and therefore for the various points of the input frame from which the energy content frame has 
been derived by the operations summarized above). For this purpose, the coefficients of the various bands 0 - 7 are 
taken for each instant of time, in other words for each point along the time dimension (in other words, the depth di- 
mension), and an interpolating polynomial which approximates the variation of these coefficients is found. This oper- 
40 ation is shown schematically in Figs. 7D and 7E. More particularly, Fig. 7D shows the coefficients of the energy content 
frame in the various sub-bands relating to track 1 and the instant Pk. These coefficients are shown aligned with each 
other in the first diagram of Fig. 7E, where PI indicates the interpolating polynomial. 

[0057] The local spectral parameter which was required is represented (as described in greater detail below) by one 
of the coefficients of this interpolating polynomial. This parameter is indicative of the variation of the spectral charac- 

45 teristics of the signal. This is because it is related, by definition, to the properties of the polynomial which approximates 
the variation of the coefficients (representing the energy content of the signal) in each band for each point of the energy 
content frame, and therefore, indirectly, of the input frame. In other words, the spectral parameter extracted from the 
interpolating parameter PI provides a qualitative indication of the variability of the spectral characteristics of the RF at 
a given point of the input frame. 

50 [0058] Since the interpolating polynomial is calculated for each coefficient (in the various decomposition bands) of 
the energy content frame, the result of this operation will be a frame consisting of the set of spectral parameters. This 
frame, referred to below as the "spectral parameter frame", is shown in Fig. 7F, where the individual values of the 
spectral parameter are indicated by a^, where i varies from 1 to n, n being the number of tracks in the frame, and j 
varies from 1 to m, m being the number of coefficients per track. 

55 [0059] The next operation consists in the classification of the values of the spectral parameter obtained in this way. 
After it had been found that a given type of tissue (neoplastic tissue for example) caused a spectral alteration in the 
RF signal, it was established thatthis spectral alteration was manifested in an alteration of the value of the local spectral 
parameter. For a given tissue (prostate gland tissue, for example) containing a neoplasm, for example, the spectral 



7 



EP 1 341 003 A2 



parameter falls into a different class of values from that containing the same spectral parameter found for healthy tissue. 
[0060] The method therefore classifies the values taken by the spectral parameter throughout the spectral parameter 
frame, in order to determine whether, in some regions of the spectral parameter frame, this parameter statistically takes 
values which fall in the (experimentally determined) class which characterizes neoplastic tissue. The result of this 

5 operation is again a frame of coefficients, which represent the frequency with which the spectral parameter falls within 
a given predetermined class of values. Fig. 7G shows this transfer schematically: a window is moved over the spectral 
parameter frame, and, for each position of the window, the number of parameters (a^) belonging to this window and 
falling within the predetermined class is found. This value becomes the coefficient fy of the output matrix or frame, 
which is referred to below as the "population frame". 

10 [0061] An operation of determining the mean of the coefficients located within a window which is moved over the 
frame is used to derive from the population frame a frame of mean values, in which each coefficient represents the 
mean value of the coefficients of the population frame within the window for each position of the window on the pop- 
ulation frame. The coefficients of this frame are used to generate an image consisting of a color scale, whose colors 
are determined by the values of the coefficients of said frame. The image obtained in this way can be superimposed 

15 on the B-mode image on a grey scale normally generated by the echograph. The colors of the color scale indicate the 
presence or absence of tumor tissue. Superimposition on the B-mode image enables the tumor tissue to be located 
within the organ under examination. 

[0062] The operations summarized above with reference to Figs. 7 and 8 will now be described individually in greater 
detail. 

20 

Application of the time-frequency transform to the sampled and digitized signal (decomposition into frequency bands 
or sub-bands) (Fig. 7C) 

[0063] The sampling and digitizing operations are not described in detail, since they are already known. However. 
25 the characteristic operations of the method according to the present invention for obtaining the mean value frame from 
the digitized frame are analyzed below. 

[0064] In order to obtain the transform frame, the digitized signal must be decomposed into frequency bands or sub- 
bands. In each band into which the signal is decomposed, there will be a matrix or frame containing the coefficients 
of thetime-frequency transform. For each frequency band, a series of coefficients equal to the number of signal samples 

30 present along each track of the digitized frame will be obtained. 

[0065] Fig. 9 shows the decomposition algorithm which provides a cascade filter which yields the coefficients of the 
transform frame. In the illustrated example, the filtering is carried out with a sequence of high-pass and low-pass filters 
and a series of decimation operations; however, the decimation operations could be omitted. The filtering carried out 
on the digitized frame is essentially equivalent to the application of the discrete wavelet packet transform. 

35 [0066] Each stage of the cascade can be made in the form of two filters which operate in parallel on the data from 
the low-pass and high-pass branches of the preceding stage. However, since each stage has the same structure, the 
simplest solution is to provide a single pair of FIR filters and to iterate the filtering, with temporary storage of the results 
of each stage so that these can then be supplied as the input for the next filtering operation, according to the diagram 
in Fig. 10. 

40 [0067] The filter coefficients can be defined by using a Daubechies 1 6 wavelet. By way of example. Fig. 6C shows 
the variation of the Daubechies 16 wavelet and Fig. 6D shows the variation of the DWPT coefficients obtained by 
filtering the digitized and sampled RF signal of Fig. 6B. 

[0068] The number of filter coefficients is generally low, being 32 for the Daubechies 1 6 filter for example, and the 
coefficients can therefore be implemented easily in a commercially available circuit, which can normally provide filters 

45 with a limited number of coefficients, not exceeding 64 for each section, with decimation for 4 and if the filter is not of 
the linear phase type, while a maximum of 128 can be obtained with two sections in cascade. The coefficients of the 
filters which provide a wavelet transform are defined recursively and are known. For example, it is possible to use the 
coefficients defined in US-A-6066098, the content of which is incorporated in the present description. 
[0069] A further advantage of the layout in Fig. 10 is that the coefficients remain the same for each iteration, since 

50 the filters do not change. 

[0070] A disadvantage of this layout is that the output data depend on those obtained from the intermediate proc- 
esses, and consequently, if, for example, the coefficients for the fourth level of decomposition are to be calculated, it 
is necessary to pass through all the intermediate levels. It would therefore be necessary to have processing times 
which increase considerably with a rise in the desired decomposition level, owing to the excessive number of filtering 

55 operations; it would also be necessary to use an additional store to store the results of the intermediate filtering oper- 
ations, thus increasing the amount of logical components present on the card. Moreover, it should be noted that the 
hardware implementation of this type of algorithm by means of an integrated digital filter causes a progressive loss of 
accuracy on each successive iteration, since the output data are truncated to a number of bits which is smaller than 
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that used by the internal arithmetic. 

[0071] To overcome these drawbacks, it is preferable to use equivalent filters to obtain each branch, starting from 
the root of the decomposition tree, in other words to parallelize the decomposition algorithm. 

[0072] This approach requires a reprocessing of the algorithm, making use of certain properties of the filtering by 
5 decimation. In particular the following equivalence can be used: 




[0073] It is thus possible to determine an equivalent filter for each step of the algorithm and to use this filter to obtain 
the coefficients of the discrete wavelet packet transform (DWPT) of any intermediate output, starting from the customary 
input data. 

15 [0074] Forexample, if the DWPT coefficients for the output of the second stage are to befound (Fig. 9), it is necessary 
to examine all the operations carried out on the signal from the input up to the output in question. To calculate the 
coefficients of the equivalent filter for the following path: 




Output data 



35 




[0075] By applying this procedure to all possible filter paths from the root, passing through the nodes, and up to the 
leaves of the decomposition tree, we obtain what is known as the parallelization of the algorithm, as shown in Fig. 11 . 
It is immediately apparent that in this configuration the output coefficients depend directly on the input data and the 
45 processing now takes place with filters which vary as a function of the frequency band to be analyzed. 

[0076] It is therefore necessary to use a number of FIR filters equal to the number of bands under examination; this 
has the advantage of eliminating the dependence of each set of coefficients of the DWP transform on the intermediate 
data. 

[0077] The major drawback, apart from the need to use a plurality of different filters, is that the length of said filters 
50 increases progressively as the frequency of the analysis decreases, and therefore it is not possible to continue beyond 
a certain level of decomposition. However, for the purposes of the present invention, it is sufficient to decompose the 
signal to the third level, in other words into eight frequency bands, and this can be achieved with a filter length well 
within the acceptable limits. 

[0078] The operation of filtering the digitized frame as described above essentially yields an encoding in bands, in 
55 other words a series of matrices or frames containing the coefficients of the discrete wavelet packet transform for each 
track of the original frame. By decomposing the digitized signal to the third level, eight bands are obtained, for each of 
which a frame of coefficients of the transform is obtained at the output of the filtering unit, the set of frames forming 
the transform frame. Each frame consists of a plurality of tracks, each comprising a series of coefficients of the discrete 
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wavelet packet transform. Fig. 12 shows the coefficients of the discrete wavelet packet transform obtained by decom- 
position to the third level (eight bands), for a single track of the digitized frame (at the input of the filter unit). 

Calculation of the energy content frame (Fig. 7D) 

5 

[0079] The next operation performed on the coefficients forming the transform frame consists in the extraction - for 
each sub-band - of a frame whose content is a function of the signal energy (the energy content frame). This operation 
is performed for each track as follows. 

[0080] For each track of the transform frame, in each frequency band into which the digitized frame has been de- 
10 composed, the energy content is calculated; in other words a certain number of consecutive samples of each track is 
examined (Fig. 13) and their absolute value is determined, after which the absolute values are added together and the 
sum is multiplied by the time interval between the samples (T wp ), which is: 



where F c = 40 MHz, the RF signal sampling frequency. The fact that T wp is eight times the sampling time is due to the 
decimation operation performed on the samples of the eight bands. 

[0081] Figure 14 shows a block diagram of the operation of calculating the signal energy by the procedure described 
20 above. In the diagram in Fig. 14, the signal energy is calculated by taking the coefficients two by two, without super- 
imposition (as in Fig. 13(a)), along the single track of the transform frame in a given decomposition band. This method 
of calculating the energy is a compromise between high resolution and low computational cost, but it is also possible 
to use a different method, for example by examining a greater number of track samples for each step and/or by su- 
perimposing the successive steps, in other words arranging for some samples to be examined more than once. Fig. 
25 13(b) shows the operating procedure in which two coefficients are examined in each step, with a superimposition of 
one coefficient between consecutive steps, while Fig. 13(c) shows the procedure in which three coefficients are ex- 
amined in each step, with a superimposition of one coefficient between consecutive steps. 

[0082] Since the operation is performed for each track of each frame in each band into which the signal has been 
decomposed, if nc is the number of samples from which the energy is calculated, there is a reduction by a factor nc in 
30 the number of samples in the output frame. With this operation, therefore, if nc= 2, the output of each band will be an 
energy content frame which has the same n umber of tracks as the transform frame, but a different number of coefficients 
per track, being equal to half, or to half minus one, of the number of transform frame coefficients, depending on whether 
the latter are even or odd respectively. 

[0083] Fig. 15 shows the flow chart for this operation. The counters i and j represent the counters of the columns 
35 and rows of the energy content frame, in other words the number of tracks (counter i) and the number of samples or 
coefficients pertrack (counter j). Col and Row represent, respectively, the number of tracks and the number of samples 
or coefficients pertrack of the energy content frame for the signal, pin and pOut represent the pointers to the transform 
frame and to the energy content frame respectively, and *pln and *pOut indicate the corresponding locations in the 
corresponding frame. 

40 

Determination of the local statistical parameter (Figs. 7D-7E) 

[0084] The next step is to extract the local statistical parameter, representing the local variation of the spectral content, 
from the energy content frame calculated in the preceding step for each band, and to generate the spectral parameter 

45 frame, which consists of the spectral parameters calculated for each point of the energy content frame. 

[0085] For this purpose, the coefficients in the various bands relating to the same instant are taken for each track of 
the energy content frame (see Fig. 7D), and an interpolating polynomial of suitable degree (Fig. 7E) is determined and 
is used to approximate the variation of these coefficients. One of the coefficients of this polynomial represents the 
required local spectral parameter. The degree of the polynomial and the coefficient are selected experimentally on the 

50 basis of an analysis of a plurality of tissues and of the characteristics of the echographic signal which these tissues 
send in response to excitation by the transducer of the echograph. In particular, experimental analysis has shown that, 
in order to achieve good characterization of the prostate tissue, it is preferable to use a fourth degree interpolating 
polynomial and to extract the third-order coefficient of this polynomial. As shown clearly in the schematic representation 
in Fig. 7D, the result of this operation will be a frame consisting of a number of coefficients equal to the number of 

55 coefficients of the energy content frame in a band, in other words a frame which has a number of rows equal to the 
number of tracks of the input RF signal and a number of columns equal to the number of columns of the energy content 
frame (the number of samples of the digitized frame reduced by a factor nc equal to the step used in the calculation 
of the energy content of the signal and by a value dependent on the decimation carried out during the filtering of the 
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signal). 

[0086] Fig. 1 6 shows schematically the operation of determining the interpolating polynomial for the coefficients in 
bands 0-7 relating to one instant of the n-th track. 

[0087] The least squares method of approximation is useful for the calculation of the values of the coefficients of the 
5 polynomial. 

[0088] This method can be used to determine the polynomial of degree m < n-1 , given a number n of experimental 
observations y Q , y 1} , y n _ 1 corresponding to the abscissae x Qi x 1} x n _^: 

10 p ( x ) =a 0 + a i x + +a m xm 

which gives the best approximation of the detected data. 

[0089] Considering the sum S of the squares of the differences between the approximate and the observed values, 

15 

n-1 

S ■- y (a 0 + a,Xj + + a m x™-y s ) 2 



20 the coefficients must be selected in such a way that S is minimal, in order to make the polynomial approximate the 
data as closely as possible. By deriving S with respect to the coefficients and equalizing the result to zero, and giver 

n-1 

and 

30 

35 we obtain the system of equations which can be written in matrix form as: 

a = Q" 1 v 

40 where a is the vector of the coefficients of the polynomial P(x). 

[0090] Fig. 17 shows the block diagram of this operation to explain how it is performed. This module is used to 
calculate the matrix Q, its inverse Q -1 and the vector v (m_ pflV), and the value of the coefficient of the desired order 
is determined, this order being the third (a 3 ) in the present case, by finding the product of the fourth row of O" 1 and the 
vector v. This value is then multiplied by 5000 and added to 5000, in order to bring the value into the range 4957.5 - 

45 5080, which is determined experimentally and the significance of which is explained below. The vector vwas calculated 
by examining the logarithm of the coefficients of the eight bands for the same instant of time and the values of the 
corresponding abscissae. 

[0091] Fig. 18 shows in detail the flow chart representing the operation required for calculating the local spectral 
parameter for all points of the energy content frame. The Row variable represents the number of samples per track of 

50 the input frame, NumBand = 8 is the number of spectral bands (band-0 ... band 7) into which the signal has been 
encoded by the time-frequency transform, TrackBand is the track number for each of the selected bands, RelCoef is 
the value identifying the row of the matrix Q _1 to be multiplied by the vector v (RelCoef = Order of Coefficient x (Degree 
of Polynomial +1)), pin and pout are the pointers to the input and output frames (the energy content frame and the 
spectral parameter frame respectively). *pln and *pOut indicate the corresponding locations within the frame. 

55 [0092] Ultimately, on completion of the operation described above, we obtain a frame or matrix (frame of spectral 
parameters) consisting of the set of third-order coefficients of the interpolating polynomial, calculated for all coefficients 
of the energy content frame in all the spectral bands. Each coefficient is a local spectral parameter, indicating the 
characteristics of the spectrum of the RF signal at each point of the frame entering the system. 
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Determination of the distribution of the spectral parameter in the reference classes (Fig. 7G) 

[0093] It was found that the spectral parameter obtained by the procedure described above took values which could 
be classified in different classes, and that it was indicative of the tissue characteristics of the examined organ. In other 

5 words, areas affected by neoplasms, for example, are characterized by values of the spectral parameter which fall in 
a given class of values. The ranges of values defining the individual classes are determined experimentally on the 
basis of a series of experimental observations. Other diffusers or reflectors present in the examined tissue (such as 
calcifications, adenomas, etc.) give rise to values of the spectral parameter which fall into a different class, in other 
words into a different range of values. 

10 [0094] The next operation to be performed is therefore the determination of the distribution of the spectral parameter 
in the relevant classes. For example, once it has been found experimentally that prostate tumor tissue causes a spectral 
alteration in the echographic signal characterized by values of the local spectral parameterfalling within a certain range 
of classes (or set of classes) of values, the examination of the prostate will be carried out by classifying the values 
taken by the spectral parameter at the various points of each echographic image (frame), and by establishing for each 

15 region of the image the frequency with which the spectral parameter falls within the range of values (in other words in 
the class or set of classes) characterizing neoplastic tissue. 

[0095] For this purpose, a window of adjustable dimensions is defined and moved over the frame of the spectral 
parameters which has been calculated as described above. The window is moved over the spectral parameter frame 
to determine (for each position of the window) how many coefficients of the spectral parameter frame, within this 

20 window, fall within the range which defines the class or set of classes characterizing neoplastic tissue. 

[0096] For example, when it has been established that the values of the spectral parameter always fall in a range 
extending from a lower limit to an upper limit, when the range between the lower limit and the upper limit has been 
divided into four hundred classes, and finally when it has been established that in the presence of a neoplasm the 
spectral parameter statistically falls within a sub-range of the total range of values (where the sub-range comprises all 

25 the classes between an upper and a lower class), the statistical investigation consists in determining how many times 
the spectral parameter, within the observation window, falls within the range defined by the following values: 



(0.30625 X (LowerClass-1)) + 4957.5 

30 

and 



(0.30625 X UpperClass) +4957.5, 

35 

where the value 0.30625 is given by 
(UpperLimit-LowerLimit)/NumberClasses 
in which 

UpperLimit=5 0 80 
40 LowerLimit=49B7 .5 

are, respectively, the upper and lower values of the range comprising all possible values of the spectral parameter, and 
NumberClasses=40 0 

is the total number of classes into which the total range between the lower limit and the upper limit has been divided. 
[0097] The algorithm which executes this calculation req u ires the use of a counter (fIFreqClass) which is incremented 

45 by unity whenever there is a comparison with a positive outcome between the value of the spectral parameter in the 
window in question and the range of values given by the reference classes. At the end of the comparison between 
each value of the spectral parameter in the window and the reference range, the fIFreqClass counter will have taken 
a value which will be the value to be entered in an output matrix or frame, in other words in the population frame. The 
counter is then reset to zero ; so that it can perform a new check on a different region of the input frame, moving the 

50 inspection window by one step. 

[0098] This operation is performed on the whole spectral parameter frame, by moving the window over it. The size 
of the window can be 3x50: in other words, it can cover a portion of the spectral parameter frame which covers three 
tracks or rows and includes, for each of the three tracks, fifty values of the spectral parameter, corresponding to fifty 
positions along the row of the frame, in other words fifty consecutive instants of time. The window is moved by one 

55 step along the row, in other words along the direction of extension of the tracks, until all the coefficients of the first three 
tracks have been covered, after which it is returned to the start of the tracks and moved by one step along the direction 
orthogonal to the extension of the tracks. 

[0099] Fig. 1 9 summarizes this process in a block diagram, while Fig. 20 shows the detailed flow chart. The variables 
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Cof and Row represent, respectively, the number of tracks and samples per track of the output frame, pin and pOut 
represent the pointers to the input and output frames, and *pln and *pOut indicate the corresponding locations within 
the frame. Extlnf and ExtSup denote the limits of the sub-range in question, in other words the range of values in which 
the spectral parameter statistically falls when tumor tissue is present. The 3x50 window is moved over the statistical 
5 parameter frame, by one step at a time along y, in other words along the direction of each track. 50 samples of 3 tracks 
are evaluated at each step. When the window has scanned all the samples of the first three tracks, it is moved one 
step along x, the variable fiFreqClass is reset to zero and a new scan is started. 

[0100] This population matrix or frame is used directly by the display module to create an image consisting of a color 
scale which can be superimposed on the B-mode image of the echo graph. To enable the color scaleto be superimposed 

10 on the B-mode image only in the location where it has diagnostic significance, it is preferable to filter the mean values 
frame with an absolute or relative threshold which can be set by the user. Thus the B-mode image will show only those 
points for which the value of the corresponding coefficients in the mean values frame is greater than the threshold 
value. This is because the tumor phenomena which are to be highlighted are only present in these regions, while the 
color scale has no significance in the other regions, and would impede the observation of the underlying greyscale B- 

15 mode image. 

[0101] As mentioned in the introduction, processing by extraction of the energy content of the signal and interpolation 
with an interpolating polynomial to extract a spectral parameter represented by one of the coefficients of said polynomial 
is only one of the possible processes included within the broader concept on which the present invention is based. 
[0102] In a different embodiment of the invention, a direct comparison can be made between two specific frequency 

20 bands, to extract a local spectral data element based on the variation of the coefficients in these two bands. 

[0103] An implementation of this type is described below with reference to Fig. 23. F1 and F2 show schematically 
two frames of the transform, obtained by encoding the digitized and sampled RF signal into frequency bands by the 
DWP transform. In practice, the frame F1 contains the coefficients of the frequency bands from 5 to 7.5 MHz, while 
the frame F2 contains the coefficients of the frequency bands from 7.5 to 10 MHz. A window is moved over the two 

25 frames F1 and F2, this window being represented by the black square in Fig. 23 superimposed on the two frames. For 
each step, in other words for each position of the window, the maximum value M1 taken by the coefficients in the 
window on the frame F1 is determined. Similarly, for each position or step of the window on the frame F2, the maximum 
value M2 taken by the coefficients lying within the window is determined. The values M1 and M2 determined for each 
position of the window can form two intermediate frames. The first intermediate frame will consist of coefficients M1 y. 

30 while the second intermediate frame will be formed from coefficients M2y. The ratios of the coefficients of the two 
intermediate frames are found, thus generating a frame of local spectral parameters, in which each coefficient is rep- 
resented by the ratio M1/M2. The frame of local spectral parameters, indicated by F in Fig. 23, will therefore consist 
of coefficients Sij, given by: 

35 Sjj = M1jj/M2jj 

[01 04] This spectral parameter frame can be processed further by comparing the values of the individual coefficients 
with minimum and maximum threshold values determined experimentally. The coefficients lying below the minimum 
40 value and above the minimum value are rejected, while the others are represented on a monitor in the form of a color 
scale as a function of the values of said coefficients, to create an image superimposed on the B-mode image. 
[01 05] It has been found that this procedure can reveal the presence of blood vessels in the organs under examination 
without making use of the Doppler effect. 

[0106] One of the possibilities offered by the wavelet coefficient is the extraction of information on the phase of the 
45 signal. This is because the wavelet transform, by its nature, can identify even small discontinuities present in the signal, 
a mode in which there are generally small phase shifts due to a distribution of diffusers which make up the structure 
of a tissue. If the characteristics of these micro-distributions can be determined, it is possible to identify a tissue and 
distinguish between healthy tissue, in which there is a regular structure represented by a distribution of diffusers which 
may be periodic for example, and pathological tissue, in which this periodicity is disrupted. The use of the wavelet 
50 transform not only makes it possible to discriminate particular periodicities of diffusers related to the structure, but also 
provides the possibility of adaptation to the tissue under investigation, since different mother forms of wavelet can be 
chosen according to the tissue. 

[0107] Therefore, by using the method according to the present invention, it is possible to envisage the extraction 
of a local spectral parameter which is a function of the phase variation of the echographic response signal for the 
55 tissues under examination in the various frequency bands into which the digitized frame has been decomposed. 

[0108] Clearly, the drawing shows only one practical embodiment of the invention, which can be varied in its forms 
and arrangements without departure from the scope of the principle of the invention. The presence of any reference 
numbers in the attached claims has the sole purpose of facilitating reading with reference to the preceding text and 
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the attached drawings, and does not limit the scope of protection of the claims. 
Claims 

5 

1 . Method of spectral analysis of a radio frequency ultrasonic signal reflected by a structure subjected to echographic 
examination, comprising the steps of: 

a) transmitting an ultrasonic excitation signal into a portion of said structure subjected to examination; 
10 b) receiving a radio frequency response signal from said structure; 

c) applying a time-frequency transform to said radio frequency response signal, thus dividing the radio fre- 
quency response signal into a plurality of frequency bands; 

d) calculating a local spectral parameter from the values of the time-frequency transform. 

15 2. Method according to Claim 1 , in which said time-frequency transform is a wavelet transform. 

3. Method according to Claim 2, in which said time-frequency transform is a discrete wavelet packet transform 
(DWPT). 

20 4. Method according to one or more of the preceding claims, in which the radio frequency response signal from said 
structure is sampled and digitized, and the sampled and digitized signal is encoded into frequency bands. 

5. Method according to one or more of the preceding claims, comprising the steps of: 

25 - selecting at least two frequency bands into which said response signal has been divided; 

extracting said local spectral parameter from the coefficients of said at least two frequency bands. 

6. Method according to one or more of the preceding claims, comprising the steps of: 

30 - acquiring an input frame, consisting of a plurality of tracks, each track consisting of the radio frequency signal 

received from the transducer after the generation of an excitation pulse; 

sampling and digitizing each track of said frame, thus generating a digitized frame consisting of a plurality of 
radio frequency tracks, each containing a set of digitized samples; 

applying said time-frequency transform to the digitized frame, thus encoding said digitized frame into a plurality 
35 of frequency bands, generating one frame of the transform in each band; 

generating from the coefficients of at least two of said frames of the transform a frame of local spectral pa- 
rameters, in which each coefficient is obtained by processing the coefficients of said at least two frames of 
the transform. 

40 7. Method according to Claim 6 ; in which the coefficients of each of said at least two frames of the transform are 
additionally processed to obtain at least two intermediate frames, and in which the coefficients of the frame of local 
spectral parameters are obtained by processing the coefficients of said intermediate frames. 

8. Method according to Claim 7, comprising the step of determining the statistical distribution of the local spectral 
45 parameters and creating a population frame, whose coefficients are indicative of the statistical distribution of said 

local parameters in one or more reference classes. 

9. Method according to Claim 7 or 8, in which an energy content frame is generated from each of said frames of the 
transform, said energy content frames constituting the intermediate frames. 

50 

10. Method according to one or more of the preceding claims, comprising the steps of: 

acquiring an input frame, consisting of a plurality of radio frequency tracks, each track consisting of the re- 
sponse signal corresponding to an excitation signal; 
55 - sampling and digitizing each track of said frame, thus generating a digitized frame consisting of a plurality of 

radio frequency tracks, each containing a set of digitized samples; 

applying said time-frequency transform to the digitized frame, thus encoding said digitized frame into a plurality 
of frequency bands, generating one frame of the transform in each band; 
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calculating the energy content of the signal for each track of the frame of the transform in each band, thus 
defining an energy content frame, said energy content frame having a set of samples representing the energy 
content of the input signal for each track and in each band; 

interpolating for each track of the energy content frame, with an interpolating polynomial above the second 
5 degree, the samples of all the bands for each instant of time, one of the coefficients of said interpolating 

polynomial representing said local spectral parameter, and generating a frame of the spectral parameters 
containing the values of said coefficient for each coefficient of the energy content frame. 

11. Method according to Claim 10, in which said interpolating polynomial is a fourth-degree polynomial. 

10 

12. Method according to Claim 10 or 11, in which said coefficient of the interpolating polynomial is the third-order 
coefficient. 

13. Method according to Claim 10. 11 or 12, in which the statistical distribution of the coefficients of the spectral pa- 
's rameter frame in one or more classes of values is determined. 

14. Echography device comprising an echography probe and means of acquiring and processing a radio frequency 
return signal from a structure subjected to echographic examination, characterized in that said acquisition and 
processing means are programmed to implement a method according to one or more of Claims 1 to 13. 

20 
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Acquire an RF signal 
frame (input frame) 



Sample and digitize the RF signal 
and generate the digitized frame 



Apply discrete wavelet package 
transform to digitized signal and decode 
it in bands to generate the frame of 
coefficients of the transform in each 
sampling band (transform frame) 



Fig. 8 



Caiculate signal* energy 
and generate energy 
content frame for each band 



Use an interpolating polynomial to approximate 
the variation of the energy content In the various 
bands for each track of the energy content frame 



r 

"Extract the 3rd-order coefficient 
of the interpolating polynomial 



Construct the spectral parameter 
frame 



Classify the spectral parameters 
and determine the frequency frame 



Mean value of the coefficients., 
of the frequency frame and 
determine the mean values frame 
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Fig. 11 
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Fig.13 
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Fig.14 



Examine nc consecutive 
samples of a track 



Examine the next 
track 




Find the absolute values 
of the samples and add 
together 



Multiply the value of the sum 
by the time step between 
samples 



y 

END 




Enter the result in the energy 
content frame 



Move in the transform 
frame by a distance 
equal to the chosen step 
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Fig. 16 
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Calculate the matrix Q 
and its inverse Q" 1 



Fig.17 



Examine the next 
track 



yes 




Identify the samples of the 
various bands for the same 
instant 



Calculate the vector v by the 
feast squares method of 
approximation 



Calculate the value of the 
coefficient of the desired 
order 



Add 1 to the calculated value of the 
coefficient, multiply by 5000 to bring 
the value into the range 4957.5-5080, 
and enter the result in the spectral 
parameter frame 



Move by one sample (row) 
in the input frame 
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^ START ^ 



Fig. 19 



Position the window 
(dimensions 500 x 3) 
on the first 50 samples 
of the first 3 tracks 



Count the number of values within 
the window lying within the fixed 
range, increment one variable for 
each comparison with positive 
outcome 



t 

Enter the value of the 
variable in the output 
frame 




Move the window by one 
step and examine new 
samples on the same tracks 



no 




Move the window by one 
track and start again from 
the initial samples 
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^ START ^ 



Position the window 
(dimensions 8x8) 
on the first eight samples 
of the first eight tracks 



Fig. 21 



Find the sum of the values 
of the input frame contained 
in a window of fixed dimensions 



Divide the sum by the 
dimensions of the window 
to find the mean 
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Move the window by one 
step and examine new 
samples on the same 
tracks 



yes 




Move the window by one 
track and start again from 
the initial samples 
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Fig.23 



Ml= maximum over the 
highlighted range 
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to band 1 
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highlighted range 



HI 



Plot relating 
to band 2 



Each element of 
the range takes the 
va/ue M7/M2 



Matrix of 
the ratios 



38 



